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We extend our recently advanced model on collisional energy loss of 
heavy quarks in a quark gluon plasma (QGP) by including radiative energy 
loss. We discuss the approach and present calculations for PbPb collisions 
at y/s = 2.76TeV. The transverse momentum spectra, RAA, and the 
elliptic flow V2 of heavy quarks have been obtained using the model of 
Kolb and Heinz for the hydrodynamical expansion of the plasma. 

1. Introduction 

There are many pieces of evidence that in ultrarelativistic heavy ion 
collisions a plasma of quarks and gluons (QGP) is formed. Such a state of 
matter, in which unbound gluons and quarks are in local thermal equilib- 
rium, has been predicted by lattice gauge calculation, the numerical solution 
of the Lagrangian of Quantum Chromodynamics (QCD) on a lattice. 

In heavy ion reactions such a state can be only created for a very short 
time (of the order of 10 _23 s). Then the system, which expands with almost 
the speed of light, forms hadrons which are finally observed in the detector. 
The problem is now to conclude from the observed hadrons on the existence 
and the properties of such a QGP. This is all but easy. It turned out that 
the multiplicity of light hadrons is very well described by statistical models. 
This means that at the end of the plasma expansion, when hadrons are 
formed, the system is in thermal equilibrium and therefore light hadrons do 
not carry information on the early stage of the expansion. 

Therefore only few observables can be used to study the interior of the 
plasma. They include the high pt hadrons which originate from jets a as 
well as the pt and vi distribution of heavy mesons which contain either a c 
or a b quark because neither jets nor heavy quarks come to an equilibrium 
with the plasma. 
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Heavy quarks are produced in hard binary initial collisions between the 
incoming protons. Their production cross sections are known from pp col- 
lisions and can as well be calculated in pQCD calculations. Therefore the 
initial transverse momentum distribution of the heavy quarks is known. 
Comparing this distribution with that measured in heavy ion collisions al- 
low to define Raa = (d& aaI dpi) / (N c da pp /dp^), where N c is the number 
of the initial binary collisions between projectile and target. The devia- 
tion of Raa from one measures the interaction of the heavy quark with the 
plasma because the hadron cross sections of heavy mesons are small. The 
heavy quark does not come to thermal equilibrium with the QGP therefore 
Raa contains the information on the interaction of the heavy quark while 
it traverses the plasma. In addition, the distribution of heavy quarks at the 
moment of their creation is isotropic in azimuthal direction, therefore the el- 
liptic flow V2 =< cos2{(j) — (f)ji) >, where (f> ((f) r) is the azimuthal angle of the 
emitted particle (reaction plane) is 0. The observed finite v% value of the ob- 
served heavy meson can only originate from interactions between light QGP 
constituents and the heavy quarks. The simultaneous description of Raa 
and V2 and their centrality dependence, presently the only observables for 
which data exist, give then the possibilities to understand the interactions 
inside the QGP. 

Unfortunately the experimental results depend not only on the elemen- 
tary interaction but also on the description of the expansion of the QGP 
PQ . Therefore the ultimate aim is to control the expansion by results on the 
light meson sector. This has not been achieved yet for the LHC and there- 
fore it is difficult to asses the influence of the expansion on the observables. 
We use here the approach from Kolb and Heinz which has reasonably well 
described the midrapidity light mesons at RHIC [2j. We adjust only the 
charged particle multiplicity to the value measured at LHC. 

The Raa of 0.2 values observed for large pt heavy mesons are much 
smaller than originally expected. Early theoretical approaches based on 
perturbative QCD (pQCD) calculation gave much larger values and it has 
been doubted, whether pQCD is the right tool to describe this interaction. 
This early calculation, however, used ad hoc assumptions on the coupling 
constant a s and the infrared regulator /i. With a standard choice \x and 
a s an artificial K factor, an overall multiplication factor of the elastic cross 
section of around 10 [3J H] had to be introduced to match the experimental 
data. 

A while ago we advanced an approach for the collisional energy loss of 
heavy quarks in the QGP [5 } 61 [7] in which a) [i has been fixed by the demand 
that more realistic calculations using the hard thermal loop approach give 
the same energy loss as our Born type pQCD calculation and b) the coupling 
constant is running and fixed by the sum rule advanced by Dokshitzer and 
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Fig. 1. (Color online) The five matrix elements which contribute to the gluon 
bremsstrahlung. 



later used by Peshier. Both these improvements increased the cross section 
especially for small momentum transfers and reduced therefore the necessary 
K factor to 2. Here we include in addition the radiative energy loss [SJ E]. 

2. Model 

Our approach extended by including radiative energy loss has been well 
described the heavy quark data at RHIC. Therefore it is worthwhile to 
calculate what we expect for LHC energies if we modify teh model only 
in a minimal way by adjusting the initial condition do dN/dy = 1600, as 
observed at RHIC. To include radiation we have to consider the following 
5 matrix elements, displayed in fig. [H which contributes to radiation. The 
commutation relation 

T b T a = T a T b - if abc T c (2.1) 

allows us to regroup the 5 matrix elements into 3 combinations, each of 
them being independently gauge invariant: 

iM h.f° = C ai(Mi + M 2 ) 

m QED = + ^ 

iM QCD = c c i(M 1 +M 3 + M 5 ). (2.2) 
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Fig. 2. (Color online) The transverse momentum distribution of Raa at midrapidity 
for different centralities and for bottom (blue) and charm (red) quarks. Top: The 
details at small pt, bottom left (right): Raa at large p t for D-mesons (B-mesons). 



h.q. (l.q.) mark the emission of the gluon from the heavy (light quark) 
line. C a , C' a and C c are the color algebra matrix elements. The matrix el- 
ements labeled as QED are the bremsstrahlung diagrams already observed 
in Quantum Electrodynamics (QED), whereas that labeled QCD is the gen- 
uine diagram of Quantum Chromodynamics (QCD). The QCD diagram is 
the main objet of interest here because it dominates the energy loss of heavy 
quarks. We evaluate the matrix elements in scalar QCD (see ref.[TU]). They 
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are given by 

■*,SQCD n ^ (jPb+PzY n r u {Pa+Pi- k) u (2 Pa - k)e v , 
lMl = CA{19) ~^mf DAP *- Pb]{ (p a - k) 2 - m 2 £ ) 

iM SQCD = Cc{ i g f D w' [p3 _ pb]D uu' [pi _ ^ [ g ^ l(pa - pi+ps - p b ) a + 

9u'aiPi ~Pa- k)fj/ + g afl > {p b -p 3 + k) v >] e a 

(P3+PbT(Pa+Pl) U 



(P3 ~Pb) 2 (Pl -Paf 



2.3) 



M3 is obtained by replacing p a — > pb and p\ — > p% in M\. Using light cone 
gauge and keeping only the leading term in ^fs we find that the square of 
the matrix element factorizes 

\M\ 2 = \M elast (s,t)\ 2 P g (m,t,k t ,x) (2.4) 

with \M e i ast (s, t)\ 2 = <? 4 ^f- being the matrix element squared for the elastic 
cross section in a coulomb-like interaction between the heavy quark and a 
light quark (gluon). P g (m,t, s, kt) describes the distribution function of the 
produced gluons. To discuss the physics we adopt the following light cone 
vectors 
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(2.5) 



The scalar product is defined as p a Pb = Va Pb ~^ Pa Pb — PatPu and q 2 = 
t ~ q 2 . In this coordinate system x is given by k + = xp~^ and represents 
the relative longitudinal momentum fraction of the gluon with respect to 
the incoming heavy quark. In this coordinate system \Msqcd\ 2 has a very 
simple form: 



\M S qcd\ 2 = g 2 D® CD A{l-x) 2 \M elast \ 2 



h h - q t 



k 2 + x 2 m 2 _ k t ) 2 + x 2 m 2 ) 

(2.6) 
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Fig. 3. (Color onlinc)Left: Raa as a function of p t for 0-20% centrality, right: 
centrality dependence of Raa ■ We compare data from the Alice collaboration [T3"] 
with our prediction. 




Fig. 4. (Color online) pt dependence of i>2- On the left hand side we compare our 
calculations for D and B meson for [30-50%] centrality with the experimental data 
shown as this conference [] , on the right hand side we display V2 for minimum bias 
separately for c-quarks and D-mcsons and b-quarks and B-mesons, respectively. 



with the color factor D QCD = C A * C q J = §. The first term in the bracket 
describes the emission from the incoming heavy quark line, the second term 
the emission from the gluon. This shows that in light cone gauge and in 
this coordinate system in leading order of ^/s the matrix element for the 
emission from the light quark do not contribute. In the case of massless 
quarks we recover the squared matrix element of Gunion and Bertsch (GB) 
of ref.pi]. 
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3. Results 

Having the matrix elements we can calculate the cross section of the 
elastic and radiative collisions of the heavy quarks with the plasma particles. 
At RHIC we found that the agreement is best when we multiply all cross 
section with a constant K factor of 0.6. A K factor of one is also compatible 
with the data but at the limits of the error bars. These cross sections are 
embedded in the plasma expansion as described in refs. El [7]. Fig. 
[2] displays the pt dependence of Raa at midrapidity for different centrality 
bins and for c and b quarks separately. Charm quarks, being lighter, suffer a 
larger energy loss than bottom quarks and are therefore pushed more toward 
low p t . Below a centrality of 40% Raa does not change substantially. At 
small momenta we see an enhancement. There the energy loss accumulates 
the heavy quarks. For large pt, shown in the bottom part of fig. El radiative 
collisions act differently than elastic collisions. If we employ only elastic 
collisions (model E, with a K factor of 2) we see an increase of Raa with 
Pt whereas for elastic and radiative collisions (with a K-factor of 0.6) Raa 
remains almost flat. If we include the Landau Pomeranchuck Migdal effect 
which suppresses radiation we would expect a moderate increase of Raa 
with increasing pt. For comparison we display as well the calculation for 
the RHIC data which matched the experimental results. 

Fig. 3 shows the comparison of our calculations with Raa Alice data 
[13] . On the left hand side we display Raa as a function of p t of [0-20%] 
centrality. The calculations follow closely the experimental data. On the 
right hand side we see Raa for mesons with pt > 6 GeV as a function 
for the centrality. Also here we see a good agreement between theory and 
experiment 

Fig. 4 show the comparison of our calculations with recent V2 Alice data 
[14] . We see that at low pt V2 for B-mesons is substantially smaller than for 
D-mesons. This is again the consequence of the smaller mass of the c-quarks 
which can more easily absorb the V2 of the light plasma particles with whom 
they collide during the expansion. We see that the prediction of our model 
(the data have been presented for the first time at this conference when 
the calculations have been already performed) agrees with the experimental 
value in between the error bars. The right hand side highlights the difference 
of V2 between b and c quarks at intermediate pt- This difference is inherent 
in the model and may therefore serve as a verification if perturbative QCD 
is the right theory to describe the data. Whereas the V2 of D-mesons is 
slightly higher than that of the c-quarks due to the hadronisation, the heavy 
B-meson has practically the same t?2 as the b-quark before hadronisation. 

In conclusions we have shown that pQCD like models which include 
a running coupling constant as well as a infrared regulator of the gluon 
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propagator in the elastic cross section which is based on hard thermal loop 
calculations reproduce the LHC data as they reproduced the RHIC data. 
The model predicts different v-i values for D- and B-mesons as well as an 
increase of Raa in central collisions with pt for pt larger than 15 GeV. The 
model can therefore be verified by future experimental data. The results 
show that collisional as well as radiative energy loss is necessary to describe 
the data. Both contribute to Raa in a comparable way. In this analysis we 
used the hydrodynamical model of Heinz and Kolb. It remains to be seen 
how other models for the expansion of the plasma change the numerical 
values of Raa and V2- Studies of different expansion scenarios as well as of 
the influence of the Landau Pomernschuk Migdal effect are under way. 
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